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Equation for anomalous diffusion in momentum space, recently obtained in [1], 
is solved for the stationary and non-stationary cases on basis of the appropriate 
probability transition function (PTF). Consideration of diffusion for heavy particles 
in a gas of the light particles can be essentially simplified due to small ratio of the 
masses of the particles. General equation for the distribution of the light particles, 
shifted in velocity space, is also derived. For the case of anomalous diffusion in 
momentum space the closed equation is formulated for the Fourier-component of 
the momentum distribution function. The effective friction and diffusion coefficients 
are found also for the shifted distribution. If the appropriate integrals are finite 
the equations derived in the paper are applicable for both cases: the PT-function 
with the long tails and the short range PT-functions in momentum space. In the 
last case the results are equivalent to the Fokker-Planck equation. Practically the 
new results of this paper are applicable to strongly non-equilibrium physical systems. 

PACS number(s): 52.27.Lw, 52.20.Hv, 05.40.-a, 05.40.Fb 

I. INTRODUCTION 

Interest in anomalous diffusion is conditioned by a large variety of applications: semicon- 
ductors, polymers, some granular systems, plasmas in specific conditions, various objects in 
biological systems, physical-chemical systems, et cetera. 

Many concrete problems of anomalous diffusion in coordinate space have been solved on 
the basis of equations with the fractional derivatives. Recently in [2] the new approach to 
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anomalous diffusion in coordinate space has been formulated. This approach simplifies the 
problem and at the same time permits to consider the more complicated kernels (PTF in 
coordinate space). In the present paper, as well as in [1], we apply the similar approach to 
solve the problem in momentum space. 

The deviation from the linear in time < r^(t) >~ t dependence of the mean square 
displacement have been experimentally observed, in particular, under essentially non- 
equilibrium conditions or for some disordered systems. The average square separation of 
a pair of particles passively moving in a turbulent fiow grows, according to Richardson's 
law, with the third power of time [3]. For diffusion typical for glasses and related com- 
plex systems [4] the observed time dependence is slower than linear. These two types of 
anomalous diffusion obviously are characterized as superdiffusion < r^{t) >~ {a > 1) 
and subdiffusion (a < 1) [5]. For a description of these two diffusion regimes a number 
of effective models and methods have been suggested. The continuous time random walk 
(CTRW) model of Scher and MontroU [6] , leading to strongly subdiffusion behavior, provides 
a basis for understanding photoconductivity in strongly disordered and glassy semiconduc- 
tors. The Levy- flight model [7], leading to superdiffusion, describes various phenomena as 
self-diffusion in micelle systems [8], reaction and transport in polymer systems [9] and is 
applicable even to the stochastic description of financial market indices [10]. For both cases 
the so-called fractional differential equations in coordinate and time spaces are applied as 
an effective approach [11]. 

However, recently a more general approach has been suggested in [2], [12], which avoid 
the fractional differentiation, reproduce the results of the standard fractional differentiation 
method, when the last one is applicable, and permit to describe the more complicated cases 
of anomalous diffusion processes. In [13] these approach has been applied also to the diffusion 
in the time-dependent external field. 

In this paper the problem of anomalous diffusion in the momentum (velocity) space will 
be considered. In spite of formal similarity, diffusion in the momentum space is very different 
physically from the coordinate space diffusion. It is clear already because the momentum 
conservation, which take place in the momentum space has no analogy in the coordinate 
space. 

Some aspects of the anomalous diffusion in the velocity space have been investigated 
for the last decade in a few papers [14-18]. On the whole, comparing with the anomalous 
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diffusion in coordinate space, the anomalous diffusion in velocity space is weakly studied. 
The consequent way to describe the anomalous diffusion in the velocity space is, according 
to our knowledge, still absent. 

In this paper the new kinetic equation for anomalous diffusion in velocity space is derived 
(see also [1]) on the basis of the appropriate expansion of PTF (in the spirit of the approach 
suggested in [2] for the diffusion in coordinate space) and some particular problems are 
investigated on this basis. 

The diffusion in velocity space for the cases of normal and anomalous behavior of the 
PT function is presented in the Section II. Starting from the argumentation based on the 
Boltzmann type of the PTF, we derive the new kinetic equation, which in fact can be applied 
to the wide class of the PTF functions. The particular cases of anomalous diffusion for hard 
spheres collisions with the specific power-type prescribed distribution function of the light 
particles is analyzed in the Section III. The universal character of anomalous diffusion in 
velocity space is absent for this case. But for the general case of the power-type PTF with 
the different powers (which are not connected in advance) the universality takes a place. 
For this case the universal limitations for the existence of anomalous diffusion are found. In 
the Section IV the Boltzmann type equation is used to consider influence of the drift of the 
light particles on the PTF function and on the opportunity for anomalous transport of the 
heavy component. 



II. DIFFUSION IN THE VELOCITY SPACE ON THE BASIS OF A 

MASTER-TYPE EQUATION 

Let us consider now the main problem formulated in the introduction, namely, diffusion 
in velocity space {V-spa.ce) on the basis of the respective master equation, which describes 
the balance of grains coming in and out the point p at the moment t. The structure of this 
equation is formally similar to the master equation Eq. (??) in the coordinate space 

^^^^ J ^qWq,P + q)/.(p + q,^)-W^(q,p)/.(p,^)}- (1) 

Of course, for coordinate space there is no conservation law, similar to that in the momen- 
tum space. The probability transition H^(p, p') describes the probability for a grain with 
momentum p' (point p') to transfer from this point p' to the point p per unit time. The 
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(2) 



momentum transferring is equal q = p' — P- Assuming in the beginning that the charac- 
teristic displacements are small one may expand Eq. (??) and arrive at the Fokker-Planck 
form of the equation for the density distribution fg{p, t) 

dfg{P:t) ^ _d_ 

dt dp, 

Mp) = J d'qqaW{q,p); B^pip) ^^J d' qq^q^W {q, p) . (3) 

The coefficients and S^/j describing the friction force and diffusion, respectively. 

Because the velocity of heavy particles is small, the p-dependence of the PTF can be 
neglected for calculation of the diffusion, which in this case is constant Bap = SaisB, where 
B is the integral 

If to neglect the p-dependence of the PTF at all we arrive to the coefficient 74^ — 
(while the diffusion coefficient is constant). This neglecting, as well known is wrong, and the 
coefficient Aa for the Fokker-Planck equation can be determined by use the argument that 
the stationary distribution function is Maxwellian. On this way we arrive to the standard 
form of the coefficient MTA^ip) — PaB, which is one of the forms of Einstein relation. For 
the systems far from equilibrium this argument is not acceptable. 

To find the coefficients in the kinetic equation, which are applicable also to slowly de- 
creasing PT functions, let us use a more general way, based on the difference of the velocities 
of the light and heavy particles. For calculation of the function A^ we have take into account 
that the function W{q, p) is scalar and depends on q,q- p,p. Expanding W{q, p) on q • p 
one arrive to the approximate representation of the functions VF(q, p) and VF(q, p + q): 

iy(q, p) ^ Wiq) + iy'(g)(q • p) + iiy"(g)(q ■ pf. (5) 

W{ci, p + q) ~ W{q) + W'{q) (q • p) + ^W"{q) (q • p)^ + q'W'{q), (6) 

where W'{q) = dW{q, q • p)/d{qp) |q.p=o and W"{q) = d'^W{q, q • p)/9(qp)2 |q.p=o- 
Then, with the necessary accuracy, A^ equals 

Aaip) = J d'qqaqpPpW'iq) = Pa J d'qqaqcW\q) = ^J d'qq^W'{q) (7) 

If for the function VF(q, p) the equality W'{q) — W{q)/2MT is fulfilled, then we arrive to 
the usual Einstein relation 

MTA^{p)=p^B (8) 
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Let us check this relation for the Boltzmann colhsions, which are described by the PT- 
function W^(q, p) = ^^(q, p) [12]: 



. ^ 271 da 



arccos (1 - —-—),« 



/5(M2+t;2-q.v//i), (9) 



where (p = Mv) and do{x^%L)j do ^ fi and fb are respectively the differential cross-section 
for scattering, the mass and distribution function for the light particles. In Eq. (Q we took 
into account the approximate equalities for the scattering of the light and heavy particles 
= (Ap)^ = p'^{l — cosd) and 9 '2:^ Xi where p' = ixu is the momentum of the light particle 
before collision. 

For the equilibrium Maxvellian distribution /° the equality W'{q) = W{q)/2MT is evi- 
dent and we arrive to the usual Fokker-Planck equation in velocity space with the constant 
diffusion D = B/M'^ and friction (3 = B/MT = DM/T coefficients, which satisfy the 
Einstein relation. 

For some non-equilibrium situations the PTE can possess a long tail. In this case we have 
derive a generalization of the Fokker-Planck equation in spirit of the above consideration for 
the coordinate case, because the diffusion and friction coefficients in the form Eqs. (ED,© 
diverge for large q if the functions have the asymptotic behavior W{q) ~ with a < s-|-2 
and (or) W'{q) ~ l/qf^ with /3<s + 2. 

Let us insert in Eq. the expansions for W (as an example we choose s = 3, the 
arbitrary s can be considered by the similar way). With necessary accuracy we find 



dt 



dqihip + q, t)[W{q) + W\q) (q ■ p) + 



iw^"(g)(q ■ pY + q'W'iq)] - /,(p, t) [W{q) + W'iq) (q ■ p) + ^W"iq){ci ■ p)']} (10) 
After the Fourier-transformation fg{r) = J j0^exp{ipr)fg{p,t) Eq. ( |TOl) reads: 

dfg{r,t) _ f „^rw/.^ „•T^//^.^ ^„ ^ 



dt 



f ~ d 

/ dci{exp(-i(qr)[W(q) -iW'(q) (q • — ) 
J or 



-lw"iq)ic^-^r] - [Wiq)-tW'iq) (q- |:) - ^H^"(9)(q " t) (H) 

We can rewrite this equation as [1]: 
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where 



A{r) 



POO 

/ dqq"^ 


sin (qr) ^ 


Jo 


qr 



W{q) 



(13) 



5„ = r«5(r); 5(r) 



— / dqc[r[exp{—i{qr) — l]W'{q) 



% I dqq'' 



cos (qr) 



sin{qr) 
qr 



W'{q) 



(14) 



Caisir) = r„r^C(r) = I dqqaq(3[exp{-i{qr) - l]W"{q) (15) 



27r 



Cir) 



dqq^ 



2sin{qr) 2cos (qr) sin{qr) 1 



q^r" 



qr 



W"{q) 



For the isotropic function /(r) = /(r) one can rewrite Eq. (1121) in the form 

dfg{r,t) w w / ^ , ^ dfg{r,t) ^d'^fg{r,t) 
^ A{r)fg{r) + B[r)r — ^ h C[r)r 



(16) 



(17) 



dt dr dr"^ 

For the case of strongly decreasing PDF the exponent under the integrals for the functions 
yl(r), B{r) and C{r) can be expanded: 



A{r) 



dqq^W{q)] B{r) 



1 



dqq'W\q); C{r) ~ 0. 



(18) 
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Practically the approximation C(r) ^ is always applicable (see below Sec. Ill) and the 
general kinetic equation f|T2l) for the Fourier-transform of the velocity distribution function 
takes the form 



(19) 



dt 

Then the simplified kinetic equation for the case of short-range on g- variable PTF (non 
equilibrium, in general case) reads 



^ Aor fg{r} + Bor- 



dt jy\- J • ^- 

where Aq = -1/6 / dqgW(g) and Bq = -1/3 / dqq^W'{q). 
The stationary solution of Eq. (|T7j) for C(r) = reads 

, A{r') 



(20) 



fgir,t) = Cexp 



- / dr 

'o r'B{r') 



Cexp 



Apr^ 
2Bo 



(21) 
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The respective normalized stationary momentum distribution equals 

Therefore in Eq. ( 12T1) the constant C = Ng. Equation (l20l) and this distribution are the 
generalization of the Fokker-Planck case for normal diffusion on non-equilibrium situation, 
when the prescribed W^(q, p) is determined, e.g., by some non-Maxwellian distribution of 
the small particles fb. To show this by other way let us make the Fourier transformation of 
([T2|) with C = and the respective A and Ba- 

dt dp^ dpa 

Therefore we arrive to the Fokker-Planck type equation with the friction coefficient (3 = —Bq 
and diffusion coefficient D = —Aq/M"^. In general these coefficients (Eq. (fT8|) ) do not satisfy 
to the Einstein relation. 

In the case of equilibrium l^-function (e.g., = /°, see above) the equality W\q) = 
W{q)/2MTb is fulfilled. Then with necessary accuracy (the second term in Eq. f[T5]) with 
W is of order /i/M and negligible in comparison with the first one) we find A{r)/rB(r) = 
Aq/Bq = MTiy. In this case the Einstein relation between the diffusion and friction coeffi- 
cients D = (3T /M exists and the standard Fokker-Planck equation is valid. 



III. THE MODELS OF ANOMALOUS DIFFUSION IN V - SPACE 

Now we can calculate the coefficients for the models of anomalous diffusion. 

At first we calculate the simple model system of the hard spheres with the different masses 
m and M ^ m, da /do = a? /A. Let us suppose that in the model under consideration the 
small particles are described by the prescribed stationary distribution fi, = nb(f)b/uQ (where 
(pb is non-dimensional distribution, uq is the characteristic velocity for the distribution of 
the small particles) and C, = {u^ + — q ■ v//i)/MQ. 

Wa{ci,p) = ^^r d^-MO- (24) 

If the distribution (pbiO = V^^ (7 > 1) possess a long-tail we get 
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where = + i;^ — q ■ v/^)/ul. 

For the case p = the value ,^0 ^0 
anomalous W = Wa 



q /Afi ul and we arrive to the expression for 



nf,a'^TT 



Cg 
727-1 ■ 



(26) 



23-27(7 _ l)/i4-27M^-27^27-l 

To determine the structure of the transport process and the kinetic equation in the 
velocity space we have find also the functions W'{q) and W"{q). 

If p 7^ to find W'{q) and W"{q) we have use the full value ^0 = (?^/4/i^ + p^/M^ - 
q • p/Mfj,)/ul and it derivatives on q ■ p at p = 0, ^0 = — l/M/iMg and = 0. Then 



-^0 



W"{q,p) 



riba TT'y 



Therefore for p = (^0 ^0) we obtain the functions 



(27) 



W'iq) ^ 



2Mfi^ulqW ' 
The function A{r) according to Eq. ([T3] 

A{r) = A-K 



W"{q) ^ 



iAfi\ 



2„2^7+l 



riba Tcy 



2M2/X%5^27+3 



POO 

/ dqq"^ 


sin (qr) ^ 


Jo 


qr 



Wiq) = AnCa / dq 



727-3 



sin{qr) 
qr 



(28) 



(29) 



Comparing the reduced equation (see below) in the velocity space with the diffusion in 
coordinate space (27 — 1 ^ a and W{q) = C/q'^^^^) we can establish that the convergence of 
the integral in the right side of Eq. (129!) (3d case) is provided if3<27 — l<5or2<7<3. 
The inequality 7 < 3 provides the convergence for small q {q ^ 0) and the inequality 7 > 2 
provides the convergence for g — 00. 

We have establish now the conditions of convergence the integrals for B{r) and C{r) 



B(r) 



An 



dqq^ 



cos (qr) — 



sin{qr) 



qr 



W\q) 



(30) 



Convergence B{r) exists for small g if 7 < 2 and for large q —>■ 00 for 7 > 1/2. 
Finally for C{r) convergence is determined by the equalities 7 < 2 for small q and 7 > 1 
for large q 



dqq^ 



2sin{qr) 2cos {qr) sin{qr) ^ 1 



W'\q) 



(31) 



Therefore to provide convergence for A, B, C for large q we have provide convergence for 
A, that means 7 > 2. To provide convergence for small q it is enough to provide convergence 
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for B and C, that means 7 < 2. Therefore for the purely power behavior of the function fb{^) 
convergence is absent. However, for existence of the anomalous diffusion in the momentum 
space in reality the convergence for small q is always provided, e.g., by finite value of v or by 
change of the small g-behavior of W{q) (compare with the examples of anomalous diffusion in 
coordinate space [2]). Therefore, in the model under consideration, the "anomalous diffusion 
in velocity space" for the power behavior of W{q), W'{q) and W"{q) on large q exists if for 
large q the asymptotic behavior of W{q —>■ 00) ~ l/q'^"'~^ with 7 > 2. At the same time the 
expansion of the exponential function in Eqs. (fT3|) - (fT6|) under the integrals, which leads to 
the Fokker-Planck type kinetic equation is invalid for the power-type kernels W{q^,p). 

Let us consider now the formal general model for which we will not connect the functions 
W{q), W'{q) and W"{q) with the concrete form of W{q,p). In this case one can suggest 
that the functions possess the independent one from another power-type g-dependence. 

As an example, this dependence can be taken as the power type for three functions 
W{q) = a/q"', W{q) = b/q^ and W"{q) = c/q'^ , where a, (3 and rj are independent and pos- 
itive. Then as follows from the consideration above the convergence of the function W exists 
if 5 > a > 3 (for asymptotically small and large q respectively). For the function W'{q) 
the convergence condition is 5 > /5 > 2 for asymptotically small and large q respectively. 
Finally for the function W"{q) the convergence condition is 7 > 77 > 5 (for asymptotically 
small and large q respectively). 

For this example the kinetic equation Eq. (fT2l) reads 

92 



= Por°-V(r, t) + r^-^P,r~ f{v, t) + r^~^ P.^r.^^fir, t). 



where 



Ana I dCC 
'0 



2-a 



sin ( 



- 1 



(32) 



(33) 



Pi = ATTb / dCqC 
'0 



2-/3 



cosC 



sin( 



(34) 



P9 



4vrc / dCC 
'0 



4-»7 



sin( cos ( sin( 1 
~^ 2^ ^ 6 



Taking into account the isotropy in r-space we can rewrite Eq. (!32l) in the form 

dfg{r, t) 



dt 



Por"-'fir) + r^-'P,—fir,t)+r'^-'P2^fir,t), 



(35) 



(36) 
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Naturally, Eqs. (132|) . (136|) can be formally rewritten in momentum (or in velocity) space 
via the fractional derivatives of various orders. Therefore, as is easy to see, for the purely 
power behavior of the functions W{q), W'{q) and W"{q) the solution with the convergent 
coefficients exists for the powers in the intervals mentioned above. The universal type of 
anomalous diffusion in velocity space in the case under consideration, therefore, exists if 
5>a>3,5>/?>2 and 7 > 77 > 5. It appears even for the cases, when the functions 
W{q), W'{q) and W"{q) have not the short-range cutting. Of cause the general description 
is also valid for the more complicated functions W, W and W" possessing the non-power 
short range parts. 

Now let us take into account the important circumstance: usually in the problem under 
consideration there is a small parameter fi/M, which can simplify description of the velocity 
diffusion. As is easy to see, e.g., on the basis of the particular cases (e.g., Eq. fl28|) ) for 
the convergent kernels of anomalous transport the term with second space derivative in 
general equations for distributions Eqs. ( fT2l) .( fT7|) is small in comparison with the term with 
the second derivative (as well as for the case of normal diffusion in velocity space). This 
smallness is of the order of the small ratio /i/M of the mass of the particles. Therefore 
for the most physically important kernels, describing the anomalous velocity diffusion the 
term with the second space derivative can be omitted and for non-stationary anisotropic 
and isotropic cases the diffusion equation respectively reads 

*M = AW/,M + B„W^ (37) 

and 

= A{r)f,{r) + B{r)r^I^ (38) 

For the case of purely power behavior W{q) = 1/g" and W{q) = 1/q^ we have, as above, 
A{r) = Pqt'^'^ and rB{r) = Pir^~^ (with inequalities 5 > a > 3 and 5 > /9 > 2). The 
stationary solution of Eq. (138!) ) (see, also (12T1) ) for the case under consideration reads 



, A{r') 1 ^ r P^r 



a-l3+2 



/f (r) = Ce.p I- j = Ce.p h ^^,^" _ ^ ^ ,, 1 (39) 

where 5>a-/3 + 2>0. 

To find the solution in isotropic non-stationary case Eq. (1381) has to be written in the 
form 

^-^^ - B{r)r^X{r^t) = A{rl (40) 
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where X{r,t) = lnfg{r,t). The general non-stationary solution of this equation can be 
written as the sum of general solution of the homogeneous equation Y[r,t) (Eq. fj40|) in 
which A{r) is taken equals zero) 



Jro r'Bir') 



(41) 



where $ is the arbitrary function and the particular solution Z(r, t) of the non- homogeneous 
equation Eq. ( HOl) : 



Z(r,t)^/f(r 



ro r'B{r') 



(42) 



Therefore 



/,(r, t) = exp [X{r, t)] = exp [Y + Z] = L(0/f (r), 



(43) 



where L{C,) is the arbitrary function of ^, which has to be found from the initial condition 
The variable ^ (r, t) equals 



ar,t)=t+ / dr'- 



t + 



+ c, 



(44) 



r'Bir') ^ ' Pi(5-/5) 

where c is the arbitrary constant, which can be omitted due to presence of the arbitrary 
function L and the values 5>a; — /3 + 2>0, 3>5 — /3>0. The general non-stationary 
solution for the case under consideration reads 



fg{r,t) = Lit + 



„5-/3 



exp 



For 



a-/3+2 



Pi{a-(3 + 2) 



Pi(5-/3) 

The unknown function L can be found from Eq. fl45l) and the initial condition fg{r,0) 
0g(r): 



(45) 



L 



,5-/3 



exp 



(46) 



Pi{a-p + 2) 

The function (f)g{r) = J d'^pexp{ipr)fg{p,t = 0) is the Fourier-component of the initial 
distribution in momentum space. By use the notation ( = r^~^ /[Pi{5 — j3)] (what means 
r(C) = [Pi (5 - /9)C]^/(^-^)) we find 



HO = (i)g[r iO] exp 



Pi{a-(3 + 2) 



(47) 
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Therefore the time-dependent solution is equal 



/,(r, t) = 0,[r(C + t)] exp { °^ ^ ^ , \ exp 



(48) 



Pi(a-/5 + 2)_ 

where we have express r(C + t) = [-Pi (5 — /3)(C + t)]"'^/*-^^'^-' as a function of r, t: 

r(C + t) = [Pi(5 - /3)(C + t)]'/('-^) = [r^-^ + Pi (5 - /?)t]^/(^-'^), (49) 

or finally 



fair, t) = <P, {[r'-^ + Pi(5 - /3) tj^/^^"^)) exp 



It is necessary to stress that for the fractional powers 1/(5-/5) and (or) {a — j3 + 2)/{5 — (3) 
the real solution exists only if Pi > 0. The limit t — oo for the solution only for the specific 
initial conditions can coincide with the stationary solution. 

For the power dependence of the functions W{q) and W'{q) the equation (1371) can be 
formally written in fractional derivatives: 

= Po/^V.(p,t) - i^i(3 + l)D^'f,ip,t) + Pip„DrV,(p,t), (51) 

where u — a — 3, ■y = P — 5 and D]^^ fgijp.t) = i J (Prexp{—ipr))rar'^fg{r,t). 

The specific case of anomalous diffusion in velocity space, which leads to the equation 
similar to one in [14,15] is derived in the Appendix on the basis of general equation f[T^ . 

IV. THE MODEL OF DIFFUSION ON BASIS OF THE BOLTZMANN 
COLLISIONS WITH DRIFT DISTRIBUTION 

Let us consider the simplest case of non-equilibrium, but stationary distribution fb, 
namely the shifted velocity distribution. 

The evident generalization of the PT-function for the case of the shifted velocity distri- 
bution of the light particles (with the drift velocity u^) the PT-function w^(q, p) (p = Mv) 
reads 

^2 



WB(q,p,Ud) = — I duu 



arccos(l - — ^),u 



(J'^Q Jq/2p, do 

Mu' + (v - Ud)2 - q • (v - Ud)//x). (52) 



13 



Again as in the Section II to find the coefficients in the kinetic equation let us use the way, 
based on the difference between the velocities of the light and heavy particles. At the same 
time the driven velocity Ud is not, generally speaking, small in comparison with the current 
characteristic velocities u and q/ ^ oi the small particles. 

For calculation of the function Aa we have take into account that in general the scalar 
function w^(q, p,Ud) has to be taken in the form iy(q, p, Ud) = W{q,p,Ud,l,^,ri) (here 
Z = (q ■ pd)) and expanded on ^ = (q ■ p) and r] = (Mud ■ p) = (pd • p)- In fact it is the 
expansion on velocity v, which is small in comparison with other characteristic velocities 
q/fi, u and Ud- As we showed above (for the case Ud = 0) to arrive to the simple and 
solvable equation for the distribution fg, taking into account smallness of v in comparison 
with the characteristic velocities and in comparison with (q ■ v//i)) we have approximate 
the function W{q,p, Ud, I, ^, rj) ~ W{q, Ud, I, ^, r]), because we are interested mainly the high 
value of q for anomalous transport. At the same time after this type of neglecting and 
expansion on ^ and (for Ud ^ 0) f] we can arrive for the special case of the kernels (e.g., 
purely power- type kernels on q, which are often considering for diffusion in coordinate space) 
to the divergence in some coefficients of the diffusion equation created by the region of a small 
q. This divergence is really absent for the realistic PT functions, which have cutting at small 
q. This cutting for small q has the physical reasons and is not related with the neglecting 
V = p/M in the approximation W{q, p, Ud) = W{q,p, Ud, I, ^, rj) ~ W{q, Ud, I, ^, rj). 

Let us expand W{q,p,Ud,l,C,,ri): 

W{q, p, Ud) = W{q, p, Ud, I, ^, v) - Wo{q, p, Ud, I) + dW/d^ 1^,^=0 ^ + dW/drj 1^,^=0 V + 

^d'W/Oe k,=o e + \d^Wldrf |^,,=o r,^ + d^W/did'q |^,,=o ^^,(53) 

where W^iq^p^Ud^l) = W{q,p,Ud,l,^,r]) \^^^=o. Then introducing Vi{q,p,Ud,l) = 
dW/di \^^r,=o,Ui{q,p,Ud,l) = dW/dr] |^,^=o, V"2(g,p, Wd, = \d'^W/di'^ |^,^=o, f/2(g,P, Wd, = 
^ff^W/drf 1^,^=0, W2{q,p,Ud, I) = d'^W/d^dr] |5,,,=o we can rewrite Eq. (13^ in the form 

W{q,p,Ud,l,i,ri) ~ Wo{q,p,Ud,l) + VL(q-p) + 
t^i(Pd ■ p) + V2{q ■ p)2 + f/2(pd ■ p)' + W2{(\ ■ p)(Pd ■ p) ^ 
Wo{q,p = 0,Ud,l) + Vi{q,p = 0,Md,/)(q- p) + Ui{q,p = 0,Md,/)(pd ■ p), (54) 

Finally we put p = in the coefficients of Eq. and omit the terms with second order 
derivatives due to existence of the small parameter n/M. It means only the terms of order 
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■s/JifM are essential in the expansion of PF(q, p, Ud). Let us calculate W{c{, P + q, Ud) taking 
into account the difference between the values of the characteristic momenta q, p + q, Ud or 
velocities. 

Then the expansion for the function W{(\, p + q, Ud) with the necessary accuracy reads 



1 92 

iy(q, p + q, Ud) ~ iy(q, p, Ud) + {qad/dpa + -g^g/j ^ ^ }W{q, p, Ud) 



(55) 



or 



WQ{q,p = 0, Ud,l) + Vi{q,p = 0,Md, /)(q ■ p) + Ui{q,p = 0,Ud,l){pd ■ p) + 

Vi{q,p = 0,Mrf,/)q^ + Ui{q,p = 0, M^, /)(q ■ Pd) (56) 

W^(q, P + q, Ud) ^ Wo{q, p = 0, Ud, I) + 
Vi (g, p = 0, Ud, I) [q ■ (p + q)] + f/i (g, p = 0, Ud, I) [pd ■ (p + q)] • (57) 

Then the kinetic equation reads 

^^^^ = j t/q{ [Wo {q,p = Q,Ud,l) + Vi{q,p = Q,UdM^-v) + Ui{q,p = Q,Ud, I) (pd ■ p) 

+Vi{q,p = 0,Md,/)q^ + Ui{q,p = 0, M^, /)(q ■ pd)] 
fg{p + q, t) - [Wo{q,p = 0,Ud,l) + Vi{q, p = 0, Ud, 0(q " P) + 

Ui{q,p = 0, Ud, /)(pd ■ p)]/g(p, t)}58) 

After the Fourier-transformation fg{r) = J j0^exp{ipr) fg{p,t) Eq. ( ITOi) reads: 

df (r, t) f d 
^ = / dq{exp{-iqr)[Wo{q,p = 0,Ud,l) -iVi{q,p = 0,Ud,l){ci- —) 

d 

-iUi{q,p = 0,Ud,l){pd ■ ^)]/9(r,i)) 

d d 
-[Wo{q,p = 0,Ud,l) -iVi{q,p = 0,Ud,l){q- ^) - iUi{q,p = 0,Ud,l){pd ■ T-)]fg{r,t))} 

or or 

(59) 

Therefore 

^^^M = Ad{r)fg{r,t) + Bd{r)^fg{r,t), (60) 

where 

Ad{r,pd) = J dq[exp{-iqr) - l]Wo{q,p = 0,Ud,l), (61) 
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Bd(r, pd) = -i J dq[exp{-i(ir) - l]{Vi{q,p = 0, m^, /)q + = 0,Ud,l)pd} = 

ri?^(r,pd) + Pd5^'(r,Pd), (62) 

In more detail we consider this equation in the separate paper. 



V. APPENDIX. ANOMALOUS VELOCITY DIFFUSION FOR THE SPECIFIC 

CASE B(R)=CONST, C(R)=0 

Let us consider now formally the specific particular case of anomalous diffusion, when 
the specific structure of the PTF (q, p) provides a rapid (let say, exponential) decrease 
of the functions W'{q) and W"{q). Therefore the exponential function under the integrals 
in the coefficients B{r) and C{r) can be expanded, that means B{r) = Bq and C(r) ~ 
respectively. At the same time the function W{q) = a/q"" has a purely power dependence 
on q. 

Then the kinetic equation Eq. f|T2l) reads 

= Por"-V.(r,t) + Sor.|-/,(r,t), (63) 

or formally in the momentum space 

= PoD^Uv.t) - Bo^^[pd,{p,t)], (64) 

where z/ = (a — 3) {2 > u > 0) and we introduced the fractional differentiation operator 
D'^f{p,t) = J drr" exp{—ipr) f {r , t) in the momentum space to compare this equation with 
the similar one in [14]. The stationary solution of Eq. (jHS]) is equal 

/,(r) = Cexp [-^^] (65) 

/g(p) = /" d^rexp{—ipr)exp [ ^ ] = — ^ [ drrsin{pr)exp [ ] (66) 

J Bo p Jo Bo 

For the case z/ = 1 we find /(p) = ng6{p), therefore C = UgCxp [Po/-Bo]/(2vr)'^. The similar 
consideration has to be used for other types of anomalous diffusion in velocity space. The 
physically important applications based on the physical models for W{(\, p) function will be 
considered separately. 
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VI. CONCLUSIONS 

In this paper the problem of anomalous diffusion in momentum (velocity) space is con- 
sequently considered. The new kinetic equation for anomalous diffusion in velocity space is 
derived without suggestion about existence of the equilibrium stationary distribution func- 
tion. Namely for the strongly non-equilibrium situations the long tails in PT-functions can 
manifest themselves. The model of anomalous diffusion in velocity space is described on the 
basis of the respective expansion of the kernel in master equation. The conditions of the 
convergence for the coefficients of the kinetic equation are found for the particular cases. The 
wide variety of the anomalous processes in velocity space exists, because even in isotropic 
case the three different coefficients in the general diffusion equation are present (one of them 
is usually negligible) . The example of the Boltzmann kernel with the prescribed distribution 
function for the light particles is studied, in particular for the hard spheres interaction. In 
general the Einstein relation for such situation is not applicable, because the stationary state 
can be far from equilibrium. For the normal diffusion the friction and diffusion coefficient are 
explicitly found for the non-equilibrium case. For equilibrium case the usual Fokker-Plank 
equation is reproduced as the particular case. 

The non-stationary and in general non-equilibrium for i — > solution is found for the 
definite initial conditions. 

The kinetic equation for the heavy particles distribution in the case of the prescribed 
distribution function for the light particles, which possess a drift velocity is derived in the 
suggested general approach. 

Acknowledgment 

The author is thankful to E. Allahyarov, W. Ebeling, M.Yu. Romanovsky and I.M. 
Sokolov for valuable discussions of some problems, reflected in this work and E. Barkai for 
the useful remark. I also express my gratitude to the Netherlands Organization for Scientific 
Research (grant NWO 047.017.2006.007) and the Russian Foundation for Basic Research for 



17 



support of my investigations on the problems of stochastic transport. 



[1] S.A. Trigger, ArXiv 0907.2793 vl, [cond-matt. stat.-mech.], 16 July 2009; in print. 
[2] S.A. Trigger, G.J.F. van Heijst and P.P.J. M. Schram, Journal of Physics, Conference Series 
11, 37 (2005). 

[3] A.S. Monin and A.M. Yaglom, 1975 Statistical Fluid Mechanics: Mechanics of Turbulence Vol 

11 (Cambridge: MIT Press) 
[4] B. Rinn, P. Mass and J.-P. Bouchaud 2000 Phys. Rev. Lett. 84, 5405 (2000). 
[5] R. Metzler, and J. Klafter, Phys. Rep. 339, 1 (2000). 
[6] H. Scher and E.W. Montroll, Phys. Rev. B12, 2455 (1975) 
[7] J. Klafter, M.F. Schlezinger and G. Zumofen Phys. Today 49, 33 (1996) 
[8] A. Ott, J.-P. Bouchaud, D. Langevin and W. Urbach, Phys. Rev. Lett. 65, 2201 (1994) 
[9] l.M. Sokolov, J. May and A. Blumen, Phys. Rev. Lett. 79, 857 (1997) 
[10] P. Gopikrishnan, V. Plerou, L.A.N. Amaral, M. Meyer and H.E. Stanley, Phys. Rev. E60, 

5305 (1999) 

[11] B.J. West, M. Bologna, and P.Grigolini, Physics of Fractal Operators, (Springer- Verlag New 
York), 2003. 

[12] S.A. Trigger, G.J.F. van Heijst and P.P.J.M. Schram, Physica A 347, 77 (2005). 
[13] S.A. Trigger, Physics Letters A 372, 8, p. 1288 (2008) 
[14] H.C. Fogedby Phys. Rev. Lett. 73, 2517 (1994) 

[15] S. Jespersen, R. Metzler, H.C. Fogedby, ArXiv 9810176 vl, [cond-matt. stat.-mech.], 15 Oc- 
tober 1998. 
[16] E. Barkai, Phys. Rev. E60, 5305 (2003) 
[17] E. Barkai, Journal of statistical physics, 115, 1537 (2004) 

[18] W. Ebeling and M. Yu. Romanovsky, Microfields, kinetic equations and fusion rates in ex- 
ploding ions clusters, Contr. PI. Physics (in print, 2009). 
[19] D. Brockmann, and l.M. Sokolov, Chemical Physics 284, 409 (2002). 



